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We model the one-dimension (ID) to three-dimension (3D) crossover in a cylindrically trapped 
Fermi gas with attractive interactions and spin-imbalance. We calculate the mean-field phase dia¬ 
gram, and study the relative stability of exotic superfluid phases as a function of interaction strength 
and temperature. For weak interactions and low density, we find ID-like behavior, which repeats 
as a function of the chemical potential as new channels open. For strong interactions, mixing of 
single-particle levels gives 3D-like behavior at all densities. Furthermore, we map the system to an 
effective ID model, finding significant density dependence of the effective ID scattering length. 
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I. INTRODUCTION 

Spin-imbalanced Fermi gases are predicted to display 
an array of exotic superconducting phases, where the or¬ 
der parameter has non-trivial structure [1-24]. Mean- 
field theories predict that these states occupy a very 
small fraction of the phase diagram in 3D, but are ubiq¬ 
uitous in ID [1-9], with the caveat that quantum fluctu¬ 
ations prevent long-range order in ID [25]. Indeed, cold- 
atom experiments in 3D [27] have found no sign of the 
exotic Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase 
[26], while experiments on ID tubes [10] found thermody¬ 
namic evidence for a fluctuating version [11-15] of FFLO, 
but were unable to measure the order parameter. One av¬ 
enue for directly observing these exotic superfluid states 
is to use highly anisotropic quasi-lD geometries where 
they should be robust [14, 28-36]. Here we solve the 
Bogoliubov-de-Gennes (BdG) equations in such a geom¬ 
etry. We find large regions of the phase diagram in which 
the FFLO finite-momentum-pairing state is stable. We 
also find a stable breached-pair (BP) state where pairs 
coexist with a Fermi surface [1, 2, 37-40]. Our analysis 
provides a much needed narrative for thinking about the 
1D-3D crossover, going beyond the existing single-band 
models [28-30, 41, 42] and studies of finite systems [33- 
35]. While we focus on cold atoms, these considerations 
are also relevant to nuclear, astrophysical [40, 43, 44], 
and condensed-matter systems [44-46]. Evidence of the 
FFLO phase has recently been found in a quasi-2D su¬ 
perconductor [46] , and there are ongoing attempts to see 
related physics in 2D atomic systems [47]. 

We consider a harmonic oscillator potential of fre¬ 
quency uj± which confines the motion of the atoms in the 
x-y plane. The atoms are free to move in the z direction, 
have mass m, and interact via s-wave collisions, charac¬ 
terized by a scattering length Ug (tuned via a Feshbach 
resonance [18, 48]). We consider the “Bardeen-Gooper- 
Schrieffer (BCS) side” of resonance where Og < 0, and 
calculate the mean-field phase diagram in the y-h plane, 
where y = (/r^ -I- and h = — ui)/2 denote re¬ 

spectively the average chemical potential and the chem¬ 
ical potential difference of the two spins. Prior work on 
this model has examined the low-density (small fi) limit. 


where the transverse motion of the atoms are confined to 
the lowest oscillator level [9-11, 28, 41], and the system 
maps onto an effective ID model [49]. Gonversely, when 
y, is large, the atoms can access many energy levels of the 
trap, and the system is locally three-dimensional. Here 
we investigate the crossover between these regimes. 

The exact ID phase diagram contains three phases 
which are fluctuating analogs of the BGS superfluid, 
the FFLO state, and a fully polarized (FP) gas [9-12]. 
Since interaction effects in ID are stronger at low den¬ 
sities, pairs are more stable at smaller /x, and the slope 
7 = dfi/dh of the line separating the BCS and FFLO 
phases has a negative slope. The analogous phase bound¬ 
ary in 3D has a positive slope, providing a convenient 
distinction between ID-like and 3D-like behavior. In 3D 
there is also a partially polarized Normal (N) state [1, 2]. 

For weak interactions and fj- < 21tw±, we find ID-like 
behavior, in that 7 < 0. The critical field jumps when¬ 
ever a new channel opens (near y ~ nhuj±), but after this 
jump we again find 7 < 0 (Fig. la). Once many channels 
are occupied we find 3D-hke behavior with 7 > 0. Each 
ID-like interval hosts a large FFLO region. In the 3D 
regime, these regions merge to form a single domain. As 
interactions are increased, the crossover to 3D-like be¬ 
havior moves to smaller /i (Fig. 2a). For very strong 
interactions near unitarity (og —>■ — 00 ), the harmonic 
oscillator levels are strongly mixed, and we always find 
7 > 0. (Fig. 2b). Regardless, we find that the FFLO 
phase occupies much of the phase diagram for all interac¬ 
tion strengths. Moreover, at sufficiently strong interac¬ 
tions, we find a BP region, nestled between the BCS and 
the FFLO phases. Such a (zero-temperature) BP phase 
is stable in 3D only for negative in the deep Bose- 
Einstein condensate (BEG) side of resonance (og > 0) 
[1-5, 40, 50]. These results clearly show that exotic su¬ 
perfluids will be observable in quasi-lD experiments. 

We study the temperature variation of the phase dia¬ 
grams (Figs. 4 and 5). The FFLO and BP phases shrink 
much faster with temperature than the BCS phase as 
they have much smaller pairing energies. For weak in¬ 
teractions, the BCS phase survives in isolated pockets, 
which disappear sequentially with temperature. The crit¬ 
ical temperatures grow with interactions, as interactions 
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favor pairing. 

In addition to directly solving the 3D BdG equations, 
we map the system to an effective ID model in the single¬ 
channel limit /i < 2hu}±. We find that the effective ID 
coupling constant giB becomes more strongly attractive 
at larger g. Our mapping reduces to that in [49] in the 
low-density limit, but has previously unexplored correc¬ 
tion terms at higher densities. These become more im¬ 
portant at stronger interactions (Eq. (9)). 

We use the Bogoliubov-de-Gennes (BdG) mean-field 
formalism. This approach does not include a Hartree 
self-energy [7]. This deficiency is typically unimportant 
for weak interactions, but becomes significant as one ap¬ 
proaches unitarity. It may also be important for study¬ 
ing the competition between phases with similar energies. 
Unfortunately the literature contains no convenient way 
to incorporate the Hartree term. The technical difficulty 
is that the bare coupling constant for contact interac¬ 
tions has an ultraviolet divergence. Renormalizing this 
divergence causes the Hartree term to identically vanish, 
and there is active debate about the significance of those 
terms [7] . At unitarity, one can circumvent this problem 
by imposing universality on the equation of state, and 
constructing a regularized energy functional [22, 35, 51]. 
However, there is no equivalent scheme at intermedi¬ 
ate interactions, as the proper set of constraints are un¬ 
known. Along with self-energy corrections, quantum fluc¬ 
tuations also become significant at stronger interactions 
[52]. Thus we do not expect our results to be accurate 
in the unitary regime. In fact, a recent experiment with 
®Li atoms, performed near unitarity, found the behavior 
of the system to be ID-like at low densities [10], whereas 
our model predicts 3D-like physics there. We believe the 
physics neglected in the BdG approach largely renormal¬ 
izes Os, and that our unitary results should agree with 
experiments for Og > 0. Lastly, we cannot rule out other 
phases not considered here, e.g., a state with deformed 
Fermi surface pairing [16], or an incoherent mixture of 
paired and unpaired fermions [17]. 

Despite these limitations, our simple model lets us 
make concrete predictions, and provides insight into the 
nature of the dimensional crossover. In particular, we 
find that the phase diagram changes dramatically with 
interaction strength (Figs. 1 and 2). These phase dia¬ 
grams, and even the equation of state, can be probed in 
experiments [1, 3, 15-17, 20, 29, 31, 32, 37, 55]. 


II. MODEL 


Our starting point is the many-body Hamiltonian 


H = 


d^r[ 


cr=t,4 


■ 9 


( 1 ) 


where denote the fermion field operators, is 

the single-particle Hamiltonian, = —1 i^V^/(2to) -|- 
{l/2)muj‘j_{x‘^ +y‘^), and g is the ‘bare’ coupling constant 
describing interactions between an f-spin and a ),-spin. 
We can relate g to Ug by the Lippmann-Schwinger equa¬ 
tion l/g = m/{4:Trh^as) — Jd^km/{8TT^h^k'^) [56]. We de¬ 
fine the pairing field A('r) = g{'ip^{r)il)^{r)), and ignore 
quadratic fluctuations, arriving at the mean-field Hamil¬ 
tonian 


^MF 


J V A*(r) 

+ Yiel^ - di) - a/g) f rfV|A(rO|^ 

n 


( 

( 2 ) 


where denote the single-particle energies. We diag¬ 
onalize by a Bogoliubov transformation, obtaining 
ijMP = + (F„ ++ (e„ - A„)] 

-(1/g) J d^r|A(r)|^. Here e„ = e^P - g., ^ represent 

the Bogoliubov quasiparticle annihilation operators, and 
the eigenvalues En (^ 0) are determined from 


\ A*{r) g-H^p)\vlr)) \v{r)J ' 


( 3 ) 


In the zero-temperature ground state, all quasiparticle 
states with a negative energy are filled, and others are 
empty, which yields a total energy 


n 

where a{x) = x for a; < 0, and 0 for a: > 0. The ground- 
state solution is found by minimizing £ as a function of 
A(r) for a given g and h. 

To simplify calculations, we take the ansatz A(r) = 
Ao exp [—(a:^ -I- 2 /^)/^^] exp (* 9 z), and minimize Eq. (4) 
with respect to Ag, and q. The exp{iqz) factor de¬ 
scribes Fulde-Ferrell (FF) pairing at wave-vector q. The 
ansatz (with q = 0) also encompasses the BCS and the 
BP phases, and when Ag = 0 includes the Normal phase. 
A Larkin-Ovchinnikov (LO) ansatz, in which exp (iqz) is 
replaced by cos (qz), produces very similar results. Based 
on prior calculations, one expects that further ansatzes, 
such as the liquid crystal phases in [24], will also give 
similar boundaries. While we label regions of the phase 
diagram as FFLO, the exact nature of the order is un¬ 
certain. 

We diagonalize Eq. (3) by expanding u{r} and v{r) 
in the single-particle states with energies lower than a 
cut-off Ec- We exactly solve this finite-dimensional low- 
energy sector, and calculate the contribution of higher- 
energy states perturbatively. We write Eq. (3) in 
the bra-ket notation, and express |z)) in terms of |it) 
to obtain (iJ"P - g)\u) + A(77"p + E - g)-^A^u) = 
E\u). The second term acts as a perturbation, yield¬ 
ing En — En = {n\A{H^P + e^P — 2g)~^A^n), where |n) 
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FIG. 1. (Color online) Zero-temperature phase diagram of 
a two-component Fermi gas in a 2D harmonic trap of fre¬ 
quency uj±. Here d±/as = —3, where as is the 3D scattering 
length, and dx = . a: Phase boundaries calcu¬ 

lated using 3D BdG equations, b: BCS critical held of the full 
model (solid curve) and of various effective ID models (dashed 
curves). Short-dashed (red): ID BdG with the mapping in 
Eq. (9), dot-dashed (blue): ID BdG with Olshanii’s mapping 
[49], long-dashed (green): Bethe ansatz with Eq. (9). 

is the corresponding single-particle state. Using com¬ 
pleteness of the single-particle states, we write this as 
En — En = /g°° which 

can be expanded in powers of £~^ using the Hadamard 
lemma. Since £„ is large, we only retain the first term, 
which is (n|AAf |n)/(2£„). Thus we rewrite Eq. (4) as 

£ = £e^-'^{n\AA^\n)/{2£n) - g~^ J d^r|A(r)p, (5) 

where £ex denotes the exact-diagonalized part, and the 
sum is over n with e^P > E^.. We take |n) = \nx,ny,k), 
where {nx,ny) labels harmonic oscillator states in the 
x-y plane, and k labels plane waves along z. Then 
+ Uy + l)hu}± + h^k"^/{2m), and (n|AA'f'|n) = 
AoC^/(47rc?5_ ^/n^riy) for large n^, Uy, where (£]_ = 

h/{mui±). Thus the energy per unit length along z is 

f = £ex+A^fV(47r)[fcc(l+/((l-/i)Ac))-^/(2as)], (6) 

where kc = {2{Ec — f{x) = -s/Ar and 

the tildes denote non-dimensionalized quantites, with en¬ 
ergies rescaled by hLu± and lengths rescaled by d±. We 
perform calculations with E^ = 10. We verified that our 
results are unchanged if Ec is made larger. Our approach 
to including high-energy modes eliminates the ultraviolet 
divergence associated with the contact interaction. It is 
similar to the approach in [7], where higher modes are 
included via a local-density approximation. Other regu¬ 
larization schemes have also been successful [57]. 

III. RESULTS OF THE FULL MODEL 

Figure la shows the phase diagram at weak interac¬ 
tions. For small h the ground state is a fully paired BCS 
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FIG. 2. (Color online) Zero-temperature phase diagram of 
the full model for a: du/us = —3/2, b: d±/as = 0. Dashed 
curves plot the BCS critical field predicted by effective ID 
models. Conventions for the curves are same as in Fig. 1. 
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FIG. 3. (Color online) Quasiparticle dispersion for a: the BP 
phase at {ji,h) = (2.2,1), b: the FFLO phase at {jl,h) = 
(2,0.8), c: the FFLO phase at {jl,h) = (2.35,1.15). Different 
curves denote different transverse modes. 


state. Increasing h drives a first-order transition to an 
FFLO or a Normal (N) region. As described earlier, in 
this weak-coupling limit, the phase boundary is reminis¬ 
cent of ID, with a structure that repeats with y as various 
channels open. The FFLO state is most stable when y is 
just above ntujj± for integer n. The length ^ over which 
A{r) falls off increases with /i. The FFLO wave-vector 
q grows with h. The FFLO-Normal and FFLO-FP tran¬ 
sitions are second-order, with the amplitude Aq — > 0 as 
the boundary is approached. 

Figure 2 shows how the phase diagram changes at 
stronger interactions. As interactions favor pairing, we 
find superfluidity over a larger area. However, the phase 
diagram becomes more 3D-like, and the relative stabili¬ 
ties of different superfluid phases change. In particular, 
we see the appearance of a stable breached-pair (BP) 
phase near unitarity. As seen in the excitation spectra 
in Fig. 3a, the BP state is a gapless superfluid with a 
uniform order-parameter (in the z direction), which con¬ 
tains both paired and unpaired modes. The unpaired 
fermions fill the sea of negative energy states. The litera¬ 
ture (mostly on isotropic systems) distinguishes between 
BP states by the topology of the Fermi sea [37, 40, 50]. 
For a given transverse quantum number, the Fermi sea 
in Fig. 3a is connected, making our state analogous to 
the “BPl” state in [37]. We do not find BP states where 
a Fermi sea is broken into disjoint momentum-intervals 
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FIG. 4. (Color online) Variation of the superfiuid regions 
with temperature for dx/cis = —3. a: FFLO region, b: BCS 
region(s). The BCS phase is stable to the left of the curve(s) 
at a given temperature. 



FIG. 5. (Color online) Variation of the superfluid regions with 
temperature for d±/as = 0. a: FFLO region, b: BCS region. 


(cf. [1-5, 19, 33, 38, 40, 50, 53, 54]). However, we do find 
FFLO states of both varieties (Fig. 3b-c). The BCS-BP 
transition, as well as the BP-FFLO transition are first- 
order, accompanied by jumps in the polarization. 

We show the phase diagrams at finite temperature 
in Figs. 4 and 5. Here we include thermal fluctua¬ 
tions at temperature T by minimizing the mean-field 
free energy F = £ — TS, where S denotes the entropy. 
This has the effect of changing the sum in Eq. (4) 
to (-l//3)I]„ln(l -h -h + h), where 

P = l/(/cBr), and En takes on both positive and neg¬ 
ative values. Such a mean-field approach ignores the 
contribution of non-condensed pairs, and overestimates 
the critical temperature [3, 4, 21, 52, 58]. However, we 
expect the qualitative features in Figs. 4 and 5 to be 
valid. In particular, we find vastly different critical tem¬ 
peratures for the FFLO and BCS phases, requiring sep¬ 
arate figures to show the behavior. This separation of 
scales is reasonable, as the pairing energy of the gapped 
BCS phase is much larger than the gapless FFLO or 
BP phases. The critical temperatures grow with the in¬ 
teraction strength since the pairing energy is increased. 
The BCS phase acquires polarization at finite T, which 
causes Aq to decrease with h, making the BCS-Normal 
transition second-order at small fi. At sufficiently high 
temperature the BP and BCS phases merge and become 


indistinguishable. The most striking feature of the weak- 
coupling phase diagram (Fig. 4) is that the BCS region 
breaks up into a series of disconnected lobes which dis¬ 
appear one by one at higher temperatures. 


IV. DERIVATION AND COMPARISON WITH 
AN EFFECTIVE ID MODEL 


To further understand this system, we take q = 0 
and map it onto an effective ID model for p, < 2. We 
project Eq. (3) into the harmonic oscillator basis, treat¬ 
ing A^^fi = (m|A|n) as a perturbation if n or m 7 ^ 0 
(where n = (jix, ny)). This yields a ID BdG equation for 
the n = 0 mode. Neglecting the influence of higher modes 
on the lowest mode yields an energy per unit length 

^ _A ^ f d k V d - 2 f dk 

^ = + £fl,+ 

m,n 



where the integrals are over all fc, and the prime on the 
sum stands for (m, n) ^ (0,0). Here £5 + = -I- AS g) 2 ^ 

and £a,+ = £n for ft 0, with = Ux+Uy + k^/2 + 1 — p. 
The first two terms in Eq. (7) separately diverge, but 
their sum is finite. This expression for £ maps to that 
of a purely ID mean-field model provided we identify 
the effective ID order parameter Aid and the coupling 
constant giu as Aid = ^00 ~ + 1 ), and 


1 

fflD 


e; 

m,n 


— lim 

ric—^'OO 


c c 


TnyiTin, 


' (P 1)2 

T ’ J2nc-2p + 3 
dk Eifi !T 


Fk 




■ ( 8 ) 


Here gm = gml{d^hwi_), and Crnn = (2/(|^ + 1))"*+” 
(F(^^D^) 2 Ei(-m, -n; /{-Km\n\) when 

m -I- n is even, and 0 otherwise, 2 F 1 being a hyperge¬ 
ometric function. The prime on the sum now stands for 
2 ^ rux EUx +my -\-ny ^ 2nc. The expression in Eq. ( 8 ) 
converges as ric The effective coupling constant (?id 
is weakly dependent on A(r), and its structure is best 
understood by taking Ag —>■ 0, ^ 1, for which 


1 

5 lD 


1 

2cia 


ah‘ 2 -d) 


2 V 2 




2-21 


^ Vj +1 ~ M 


tan 


/ M -1 

j+l-A’ 


(9) 


where C denotes the Hurwitz zeta function, and 0 is 
the unit step function. At /i = 1, I/ 51 D = l/(2as) -|- 
C(l/2)/(2-\/2), which is Olshanii’s two-particle result [49]. 
As p grows, giu decreases, approaching —00 as /i —)■ 2 . 
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This divergence is unphysical, and signals a breakdown 
of the mapping to ID when more channels open. 


In Fig. lb we evaluate the validity of this mapping 
by plotting the critical field of the BCS phase, he, from 
the effective ID model. It closely follows the critical field 
obtained from the full model for nearly all /I < 2. We also 
plot he using Olshanii’s mapping [49], which agrees with 
the full model at small jl, but becomes less accurate as fl 
increases. Further, we show the prediction of the Bethe 
Ansatz with the mapping in Eq. (9), which illustrates the 
difference between an exact and a mean-field analysis in 
ID [9]. The mapping to ID becomes less accurate at 
stronger interactions due to mixing of the trap levels, as 
seen in Fig. 2. 


V. OUTLOOK 

Achieving the temperatures required to directly ob¬ 
serve the FFLO state at weak coupling is extremely chal¬ 
lenging. The numbers near unitarity are more promis¬ 
ing, but the accuracy of our mean-field theory is ques¬ 
tionable there. The ID thermodynamic measurements 
are promising [10]. Time-dependent BdG calculations 
suggest the FFLO domain walls will be observable in 
time-of-flight expansion of ID gases [31]. This signature 
should be even more robust in the geometries we have 
been studying. There are also interesting connections to 
experiments on domain walls in highly elongated traps 
[59]. It is likely that these various research directions 
will converge in the near future. 
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